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TECHNICAL NOTE NO. 859 



ON THE DISTRIBUTION OE STRESS IN A THIN PLATE 
ELASTICALLY SUPPORTED ALONG- TWO EDGES AT 
LOADS BEYOND THE STABILITY LIMIT 
By Louis G. Dunn 



SUMMARY 



This report presents the results of an investigation 
on the distribution of stress in a plate elastically sup- 
ported along two edges and subjected to compressive loadB 
in a direction parallel to the elastically supported edges 

In section I methods are developed for calculating: 

1. The compression stresses at the median surface in 

directions parallel and perpendicular to the 
applied load. 

2. Bending stresses ■ - 

3. Combined bending and compression stresses 

4. The effective width of sheet acting with, each 

stringer 

Section II presents the results of a limited number 
of tests on: 

1. The wave form of the buckled plate with increas- 

ing stringer stress 

2. Compression tests on 24S-T dural sheet and 24S-T 

alclad sheet 

A comparison between the assumed wave form and the 
experimentally determined wave forms given in reference 4 
is indicated in section I. The calculated strains have 
been compared -with the measured strains of reference 4. 
The discrepancies between theory and experiment are dis- 
cussed in sections I and II. 
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INTRODUCTION 



A considerable amount of theoretical and experimental 
work has been done on- the stress distribution in thin 
plates at loads beyond the stability limit, as is evidenced 
by the references included at the end of this report. The 
two important -factors which enter into the theoretical cal- 
culations are the amplitude and the wave form of the buck- 
led plate as a function of the ratio of the edge stress to 
the buckling stress. Several approximate expressions for 
the amplitude have been derived, (See references 1, 2, 
and 5.) Calculations of the exact wave form after buckling 
are extremely complicated even for simple cases; for exam- 
ple, a circular plate with a simply supported edge and sub- 
jected to a uniform compression load in the plane pf the 
plate (reference 6). For this reason, in the case of rec- 
tangular plates, it is customary to assume that the wave 
form after buckling remains constant and equal fro that' at 
buckling. 

The majority of invest igalro rs have assumed a wave 
form which satisfies the boundary conditions of simple 
support along all four edges. This type of edge support 
is very seldom realized in airplane construction since the 
edges parallel to the applied load are usually, support eu 
by some form of stiffener; that is, the edges are olasti- 
cally supported. 

For the practical usage of thin plates the distribu- 
tion of irhe median fiber stress and the bending stress 
is of importance. The first is important in the effec- 
tive width calculations and the latter in the calculations 
of the maximum stresses for permanent deformations. Since 
these stresses are functions of thre—plat-e "curvature.it is 
desirable to have a reasonably close approximation of the 
correct wave form. The condition of simple support along 
the edges requires zero bending moment at the edges; 
however in the actual case the bending moments along the 
elastically supported edges are different from zero. 
Also, the maximum compression stress at the median surface 
occurs at these edges, so that fo-r the effects of combined 
bending and compression stresses it is important that the 
elastic support along the edges be considered. 

In the present report a theoretical calculation has 
been carried out in which the influence of the elastic 
support has been taken into account. No attempt was made 
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to calculate the exact wave form after buckling; it was 
assumed that it remains constant and equal to that at 
buckling. The deviation of the experimentally obtained 
•wave form from the assumed one is shown in section II to 
be quite small at relatively, high stringer stresses. 
Since the amplitude also checks the experimental values 
closely, it seems reasonable to assume that the calculated 
stresses will be in good agreement with the actual values. 

This investigation, conducted at the California Insti- 
tute of Technology, was sponsored by, and conducted with 
financial assistance from the National Advisory Committee 
for Aeronautics. 

The author wishes to acknowledge his appreciation for 
the advice given by Dr. Theodor von JCarman during the prep- 
aration of this paper and to thank Dr. Hsue-Shen Tsien for 
his willing cooperat io n , and Mr. Alfred Slater for his 
assistance. 

I- THEOEETICaL DISCUSSION Of THE FEOBLEM 

lor a plate subjected to a uniform compression load 
•along two edges and elastically supported along the edges 
parallel to the load (fig. lb), the influence of the elas- 
tic support on the wave form at buckling, and on the buck- 
ling stress has been considered in reference 5. Eor sym- 
metrical elastic supports the wave form at buckling is 
given by • • - 

w = [a cosh ay + B cos py] sin {fx (l) 

where: 




a c buckling stress, pounds per square inch 
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v Foisson's ratio 

m number of half waves 

t plate thickness, inches. 

a plate length,, inches 

b plate width, inches 

If it is. again assumed that the wave form after buckling 
remains the same as at buckling, the distribution of com- 
pression and bending stresses can be calculated. Since 
the boundary conditions at the edges are given in terms 
of displacement of the median surface of the plate, it is 
convenient to deal with displacements rather than stress. 
For the plate— the following assumptions are made: 

(a) The plate under consideration is one of the panels 
of a continuous plate attached to stiffening elements of 
equal flexural and torsional rigidity, spaced a distance b 
apart. (See figure la.) 

(b) The length of the plate is large in comparison to 
its width b, so that its behavior-. is not influenced by- 
length., ........ 

Then the boundary conditions along the elastically 
supported edges are: 

u + e s t x = 0 at y = ±£ (2) 

v = 0 at y = ±SL • (3) 
3 

w = 0 at y = ±\ (4) 
and in the x direction ~ ■ 



w = 0 at 



x = 

m 



u = 



£ . a at 
st m 



x = a. 



(5) 
(5i) 



where u and v are the displacements of the median 
plane of the plate, in directions parallel and normal to_ 
the load respectively, and e gt - Is the compressive strain 
in the stiffeners. In applying expressions (l) through 
(5') to describe the deformations of the buckled plate, the 
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conditions of zero curvature and constant displacement at 
the nodal points x = a/m of the buckle are satisfied. 
The conditions along the edge y = b/2 depend upon the 
torsional and flezural rigidity of the stiffening element 
as well as the methoff'of attachment; The influence of 
the torsional rigidity is taken into account to some ex- 
tent by the assumed wave form. Since the plate under 
consideration is a panel from a continuous plate, the v , 
component is zero at the edges y = ±b/2 by virtue of 
the continuity of the plate in the y direction. This 
condition is probably closely realized in actual structures 
such as wings where the sheet is continuous around the wing 
or where the sheet is riveted to heavy fore and aft spars. 

General expressions for the displacements u and v 
can be determined from the equations of equilibrium 



— — + — o 

3 x 3y •' 



3y + 3x 



+ -#-<>■■ (6) 



and of strain 

e x * 



3X 2 V 3x ' E V 7 J 



3u 



H. + 2J. + 3JL-3JL _ 2 (1 + v) T _ (v) 



xy 3y 3x 3x 3y. u xy 

By substituting the relations of equations (7) into equa- 
tions (6) the following two equations connecting u, v 
and w are obtained: 

cPu 



3x 3 



3y 3 1 - v 3x\3x 3y J 1 - v 3x(_V ax / 
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afz + air + i_±_v aa . ax V. ■ a_|7 a^V 

. 3y 3 • 3x 2 l-.udyVSx dy. y 1 - » 3j [\ 3* '/ 

\*yJ J + 3 x a 8y 3 X S x By ° (9) 

Substituting expression (l) for w, equations (8) 
and (9) "become 

a 3 u . afu . i + v a /an a-A . . 

ar= + a y - + rrr i"x vgi + gj; = - sin 2 ^ (k i' cosh z ^ 

+ k s » cos 20y + k 3 ' sinh ay sin |3y + k 4 » . cosh ay cos Py + k s ) 

(10) 

+ ka cosh ay ain fsy + k 3 Biih ay : co's £y + k 4 sin 2fiy] 
- I k, sinh 2 ay. + kj cosh ay sin gy + k, sinh ay cos ?y 

+ ^ sin spy] (11) 

) k, . ^g.) v 
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By inspection of equations (10) and (ll) it is apparent 
that the particular integrals can be written with undeter- 
mined coefficients in the form: 

u = sin 20x[li 1 cosh 2ay + M 3 cos 2py + U s sinh ay sin py 

+ M 4 cosh ay cos. Py + M 5 ] (13) 

v = cos 20z[lT x sinh 2ay + JT 8 sin 2py + N 3 cosh ay- sin py- 
+ U 4 sinh ay cos py] + 5 5 sinh 2ay + ff s sin 2py 

+ N 7 cosh ay sin py + H s sinh ay cos Py (13) 

The coefficients . M and T$ are determined by direct sub- 
stitution in equations (10) and (11), and are as follows: 

K =~£!L M =-r&A H _ 2AB0 g qg(l-w) 
i 16 « 16 s (a a +pP) a 

M - ~2&af(l+v) ■ K -. _ -A s (ff* + Vp 3 ) - 3?fr a -pgP ) 

* " (a 2 + P 3 ) 3 5 " ' 16 |Z5 " 

H..^ S 3 ■ , ^0 4 Ma a + 0 S )(p%v0 8 )-3 - 

16 16 "» ( a 3 + p 3 ) 3 

+ (p* ■.0")Ca B -P^)]' g,,2gfa!±Zg)- 

(a 3 +P 3 ) 3 16 a 



V = -B 8 (^if - P S ) ^ _ -ABP(UjS 2 - a*) 

■16 p 7 k(f + P 3 ) 

8 . k(0 s + p 3 ) 

The particular .solution of the differential equations 
(10) and (ll) gives the primary displacements u and V 
due to bending. The complementary solution of the homoge- 
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neons differential equations (that is, the equations ob- 
tained by equating the right-hand side of equations (10) 
and (ll) to zero) gives the induced displacements due to 
bending . 

In view of the boundary condition (5 1 ) the appropri- 
ate complementary solution can be shown to be of the form 



ii = [(Cj' + C s ') cosh 20y 



+ C 2 » 20y sinh 20yJ sin 20 x - 



e .x 
st 



v =j C i - a ' J sinh 20y 

L\ 1 + v 2 i / -j 

-Cg 1 20y cosh 20y cos 20x + e'y.y 

The general expressions for the displacements u and v 
can now be written as 

u = ( 0 1 + 0 1 ) cosh 20y + 0 • 20y B inh 20y 
J_ l a • a 

+ M. cosh 2ay + M a cos 2py + M„ sinh ay sin By 
+ M 4 cosh ay cos' 6y.+ 



v = 



sin 20 x - e s tx (14) 
+ IT sinh 2ay + N sin 26y + N cosh ay sin By 

X 3 3 

+ N 4 sinh ay cos By cos 20x + 1^ sinh 2ay + k sin 2By 

+ N cosh ay sin By + IT sinh ay cos By + & ' v y (15) 

Since the- boundary conditions indicated by expres- 
sions (2) and (3) must be satisfied for all values of x,' 
it is necessary that the terms of expressions (14) and 
(15) included in the square brackets, vanish when y = ±b/2-; 
In addition, the boundary . condition (3) requites that .the 
remaining terms of expression (15) vanish when y = ±b/2. 
Prom these conditions the constants C 4 ' and C s 1 , and 
ey « can be determined, their respective values being 
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o . = _L_ \ Jj^J^ll±2^1llll^ (cosh fa 

1 cosh 0b I- ~ v , sinh 0 b cosh 0b - 0 b 

+ 0b sinh 0b) - B 1 - 

?i sinh 0b + F a cosh 0b 
0 a \ = _ ______ — _ . — _ _ 

sinh 0b cosh 0b - 0b 



3 3-1? 

1 + V 



e'y = ■_ | w sinh ab + S sin 0b + 2T_ cosh siri |_ 
" d L 6 ' _> _> 

+ 2J sinh SJ_. cos Si 
v s 2 2 

where - 1 

P- = U 1 cosh ab + M a cos 0b + M 3 sinh sin fL_ 

+ M cosh a — • cos ^i + M 
4 2 2 5 

I = 23" sinh ab + _T sin 6b 4- ff cosh — sin |- 

2 1 S 3 <3 c 

: + _f . sinh- E_ cob Bi 
4 2 2 

The constants appearing in the expressions for u 
and v are now all evaluated in terms of the quantities, 
A, B, a, 8, and 0. Therefore, the stress at any 
point on the median surface of the plate in the direc- 
tions x and y can be calculated in terms of these 
quantities from the following "expres sions : 



1— !• |> + U£ yJ " — + u By" 



+ 1 /9w\ B , 
.2 V^V 
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where 3S_ is Young's modulus for sheet 



G E r - - - - - ' ^ s 



l-v 2 y i_u 2 [3y + Bx + 2 



Substituting in equations (16) and (17) the values of u, 
v, and w given by expressions (l), (14), and (15) gives 

e r 

0 x = ir^a"[ (5, 4 ! + ^'^ cos ^^e 1 + *V cosS <# x 

+ F a sin 2 jZSx - e st + otf ' y ] (18) 

a y = 5~Jr[ <v *4' + u V> 008 2 ^ x + J- V + U V c0 < 0 X 

+ ij 8 « sin 2 S2fx + ©» - ve I (i 9 -) 
where: " 3 -I 

T^t = ^( Ci « + C 2 f ) cosh 20y + c s • 20y sinh 20y 

+ M 1 cosh 2ay + M 2 cos 2Sy + M s sinh ay sin By 

+ M 4 cosh ay cos By + M 5 ^0 

V = Cs ' ~ Cl 0 ^ cosh 2 ^ y ** V 8inn 2 ^ 

+ 20 cosh 20y) + 2a cosh 2ay + B s 2(5 cos 2By 
+ U (a sinh ay sin By + B cosh ay cos By) 

+ 3T 4 (a cosh ay cos' By - 6 sinh ay sin By) 



2a cosh 2ay + N 6 26 cos 2By + N ? (a sinh ay sin By 
+ 6 cosh ay cos By) + N 8 (a cosh ay cos By - 6 sinh ay sin $y) ] 

I ' = g- (A cosh ay + B cos By) 2 0° 
ffg 1 = |- v Ua sinh ay - BS sin By) 3 
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Before any numerical calculations of the stresses O x 
and q„ can "be carried out it is necessary to first eval- 
uate the quantities A, B, a, {J, and 0. The constants 
A and B can he evaluated from the "boundary condition 
of expression (4) and the condition that at y = 0 and 

x = it-, the deflection w is equal to the maximum ampli- 
' 2m 

tude w Q . Jrom these conditions the following values for 
A and B are obtained: 



tco S *~- <-} 
— as — — ^ Uo) 
cos - cosh Ob J 

2 2 



cosh ab 1 

* =-o ! ^~~ZZ < 21 > 

L cos £5. - cosh J 
2 2 

To evaluate the maximum amplitude, w 0 , another 
■boundary condition is necessary. At the ptjint y = 0 
and x = a/m an e mpirical expression for the stress at' 

~ ' /a \ l/s 

the median fiber of • the form o x = O cr ( -^ r - } appears 

to be in fair agreement with the experimental evidence of 
reference 4. The comparison is shown in'figure lc. 

In the expressions for a, £ , and the only un- 

known quantity for a given panel is the ..half-wave length A. 
It is therefore desirable to determine a general expres- 
sion for X. When dealing with, a plate of large a/b 
ratio, the wave length of a buckle will not be influenced 
by the plate length, but only by the plate width and the 
elastic support along the edge's y = ±b/2. Hence, for 
large a/b ratios, the ratio of the half-wave length to 
the plate width can be calculated.- as a function of the 
parameter p,, where ■ . - 

bD _ Et g b 

C ~~ 12(l-v 3 )0 



The method of calculating jj, is given in reference 5. 

The critical buckling stress of the sheet between 
stringers was shown in reference 5 to be given by the tran- 
scendental equation — 
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- e 3 tan I ye\J/-9 S + Je^+Q 3 tanh| 9* + ^« 0 



where: „ q t s t ssi. 3 ^ 

, ^ -Si-— = 1?(1 ~ gcr 

9 = TLk 

For large values of the ratio a/b the quantity \|/ 

■ is a constant for' a given value of \x and the value of 
~b/\ is given by the- expression ' 

d9 



(22) 



or 



gi—^-Mj tan i/eM/ I 9 2 + 1 (9x|/ - 29 s ) sec a 1 A|/ - 9 s 
_ 9 2 2 v 2 2 



+ jj=~ tanh -i/9M/ + 9 a + I (Q* + 29 s ) sech s l/o^ + 6 s 
V 9\|/ + 9 s '•' 2 3 2 

- ^ = 0 (23) 

Prom equations (22) and (23), the value of b/^ or ;\/b 
and the values of \|/ as functions of (j, ■ can be determined 
and are shown in figure 2. It' should be emphasized that 
these values of \|/ are the 'minimum values corresponding 
to given values of [i. 

Now observing that, 

ccb = V B 5 + 6\|/, (3d = - 9 s , . 0"b = TL)H 

it is seen that the coefficients .M and N appearing in 
the expressions for ,u and v can be determined as func- 
tions of jjl if the coefficients are multiplied by b/w 0 . 
For example 
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M, b 



16 

&b 
16 



16 



cos §1 
2 



cos 

§b - cosh 
2 2 
cos Ik -, s 

2 



-I 8 

Ob J 



-cos 65. - cosh J 
2 2 
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(2 4a) 



(24b) 



and since $b, ab and 6b are all known for a given value 



of 



the quantity 



Mi b 
5" 



can he determined as a function 



of \x. All coefficients have been evaluated in this form 
and are shown in figures 3 to 6. The constants 0 x » and 
C s ' can he evaluated in the same manner. 

In order to simplify the necessary calculations it is 
desirable to write the functions 3? ' , 3? ' and so forth, 

4 5 

in such a form, that they can be evaluated as functions of 
the parameter It can be seen that if these functions 

ar-e multiplied by 
be written in the form 



.2 



the expressions (18) and (19) can 



a* = 



(] v s ) b s I + °° 8 + 76 + * 7 cobS ^ z: 

+ sin 3 {fz + fe'y - 5^-5* 

■ ■ J 1 -V* 

[o>3? 4 + i J B ) cos 2{2$x + 1 P s 
3 ^ + i ^ sin 8 ' (fx + eyj - 



E_ w„ s 
so 



>(35) 



+ u]? 7 cos 



VE « e st 
1 v' c 



where 



e = 

y 



b S e- 



'JL 



* 4 - 



b 2 r 4 1 



w. 



3T. 



and so forth. 



The function's 3F 4 , P 5 , 3? 6 , ]T 7 , and Jg can now be eval- 
uated for various values of y as functions of the param- 
eter (j, and are shown in figures 8 to 12 inclusive. The 
constant ey is independent of y and is shown in figure 
7 as a function of |i. 
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The amplitude' w can now be evaluated from the con- 

. . . 0 - , vi/a 

dition that at y = 0 and _ =-. a/.m the stress a- = O c -l^~j • 

Writing for € flt its equivalent £_£ and noting that B 

st w 

in equation (25) applies to the sheet, the value of 

u 

is obtained as follows: 

= -O (2JBS.Y1 ._______- jY + 7. + + + ^yl - r^- 

^ Cr ^ st / (l-^)b 2 L 4 5 6 7 (i- v a)B 8t 

or (26) 

__. _ /X] /___._ m _ v c \ Oct /_5__ 

b y _ st - Eg ^a Bt 

where: .• r ^> k a 



C st stringer stress, pounds per square inch 

E Young 1 s modulus for the stringer, pounds 'per -'square 
st inch 



E_ Young's modulus for the sheet, pounds .per square 'inch 
= |_ + * 8 + » e + I 7 + »e y ] 



s 



The expression for F 9 can be evaluated as a function of 
the parameter ' p, and is shown in figure 13. Substituting 
equations (25) and (26) gives 



c r 

- ■ t&ti % - )[<*. + * > - *« + 

1 E a 

+ 7 7 cos 2 0x + P 8 sin" 0x + 1>«^.J - — 



a 



____ (l-O _ st 

'_*_ _ £i__!____.____i N \r (us . 4 + i J)oos 2 ^ 

E e+ . E a / - V 5 



st — (27) 



(i-v 2 ) ir s Vjs 8t 



+ 1 ff 6 + . cos a 0- + 1 Fg gin s 0_ + 6 ] _ J_Eg a_t_ (28) 

J (l- y2 XE st 
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The actual load carried by the plate for a given stringer 
stress and also the effective width can now he calculated. 
The total load carried by the plate is given by 

b/s 

P = 2t. ./ Cj.- dy 



= 2* J Ox 



'f °x & * 



and the effective width is by definition 

b/s 

2 w Q t a e * = 2t 

o 

or the ratio of the effective width to the plate width is 

b/s . 

C x dy (39) 



w 6 1 r_ 



0 

+ J ' • sin 3 $x - LslJl + u5 y 
where: 13 * ri 



Substituting for a x the expression (18) and inte- 
grating gives . • 

b/a 

J dy = 3-.^'r io i co s 20x4- + F 18 » cos 8 ?lx 

i , 10 » = (l-v)^ 1 + ^—^sinh 0b + 0bC 2 '(l~v) cosh 0 b 
+( Hal + UN ^ sinh ab + C ¥5|4 + ui ") Bin p D + I s , a f* a ( gM . 
+ 6M 4 ) + VU 3 J cosh ^ sin + j^-~£_ ( a « 4 - pi£ 3 ) 

+ VH 4 j sinh cos 6& +. 0b M 5 



16 NACa Technical Note Ho. 859 



r 



3 V [ IT 5 Slnh ^ + N 6 Sil1 P b 



F I 

12 



+ N„ co ah a In + IT sinh GLk co s §£, 
2 2 8 2 2. 

a 2 "TpV a sinh 2~ C0B 2~ + p c08h T sin 2 y_ 

v [" h . a e 32. A a rr B 2 8 

= 2 I j (B ^ - A a ) + ~ a sinh ah ~ ^ sin 0b 

If the expressions for F 1 , 1? « , and so forth, are 



1 » 

13 



10 ' -11 

multiplied hy — the ahove equation can be written in 
w o 

the form 
fc/a 

E w s n a 
a dy = -JL 2— I P cos 20x + F, _ + F,_ cos 0x 

(l-V s )b 1- 11 1 



/ 



+ F sin S 0x +--22ll - E§_i§l 
13 2 J 2(l-v 3 ) 



b 2 _ -h - 

where e v = — g- e , F, _ = — ~ F ' , and bo f o rth. 
* w 0 y 11 w 0 s 1X 



Siibst it ut ing the ahove expression in equation (29) (noting 

>e. 



that a s ^ is negative and F i:i =» - -~^-) and the value 



for -°- •from equation (26), gives 

= r /oeA /a i i_ r cos 20x 



(l-V a )E st v* at , 

0x + j lB sin 8 0.] + _ Jl-L- (30) 



+ F ia cos 



r 
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As indicated by equat ion (27) , the distribution of 
the stress across the plate also varies in a lengthwise 
direction. The jfcotal load carried by the plate. is given 
by the stress distribution at the nodal lines. 7or exam- 
ple, if a stiffened plate is tested in compression, the 
total load on the plate as measured in the test, is that 
which exists at the head of the test machine, namely, at 
a nodal line. The variation of the transverse stress 
distribut ion*between nodal lines is balanced by induced 
shear stresses at the stringer, and does not influence 
the total load carried by the plate. However, the com- 
pression stress in the stringer will be influenced to 
some extent by the lengthwise variation of the plate 
stress; that is, the total load in any transverse cross 
section of a stiffened panel being constant, any varia- 
tion in the sheet load must be transferred to- the stringer. 
This means that the compression' stress in the stringer 
will be influenced by the lengthwise variation of the . 
plate stress. The extent to which the stringer stress 
will vary between nodal points of the plate, will depend 
on the ratio of the stringer cro ss- section area to the 
plate cro s s- sect ion area.' 

The transverse stress distributions at the nodal lines 
and at the maximum amplitude positions have been calculated 
in terms of the ratio ax/<7cr for p, = 0 and p, = o?, and 
are shown in figure 14. Although at the nodal. points a 
considerable difference in stress intensity occurs at cer- 
tain values, of y, in going from p = 0 to p, = «, this 
difference occurs only over a narrow region so that the 
total load carried by the plate is not greatly affected 
by the elastic support along the edges.. At a stringer 
stress of 10 times the plate buckling stress, the simply 
supported plate carries 19.5 percent less load than the 
one with clamped edges. - 

The distribution of the a~ stresses at x = 

* ... . _ 2m 

is influenced to some extent, by the induced transverse 

stresses. As shown in figure 14, at the center of the 
buckle the stress changes to tension for the edge condi- 
tion p, —' 0; whereas, for p .= =° there is no appreciable 
change. However, in the latter case the edge. stress is. 
considerably reduced; whereas, the edge stress for p = O 
changes only slightly. This transverse stress distribu- 
tion is primarily of interest in the calculations of the 
maximum combined stresses, for the maximum bending stresses 
occur along the line x = S±-. 

2m 

As previously pointed out the distribution of stress 
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at the -nodal lines gives the actual load carried by the 
plate and any variation of the stress distribut ion be- 
tween the nodal lines is balanced by shearing stresses 
at the stringer. Therefore, the effective width of sheet 
acting with each stringer would be given by equation (30), 
when 0x = 1, 2, 4, and so forth, or 

w e 1 E s f E s /a C rV /S *i *14 t*»\ 

~ - 2T1- z*r ■ st - Lr^r^t \a~J J ir ( } 

where F, „ = F ^ + F ^ 

X 4 10 IS 

The values of F 14 , for various values of jj, , are shown 

in figure 15. If the calculated ef f ect ive- widt h is to be 
compared with experimental results, further consideration 
must be given to the method in which the experimental re- 
sults are obtained. From the previous discussion it is 
obvious that the compression stress in the stringer is not 
constant along its length but is a minimum air the nodal 
lines of the plate, and a maximum at the maximum amplitude 
of the plate. To obtain the experimental values of effec- 
tive width. It is only necessary to measure the stringer 
stress corresponding to a given total panel load. The 
load carried by the plate alohe can then be determined, 
and by the definition of the e-ffective width, its valuo 
becomes immediately calculable. However, if the stringer 
stress is not constant along its length, the effective 
width thus determined would be a function of 'the position 
at which the stringer stress is measured. The experimen- 
tal results given in reference 5 were obtained by averag- 
ing the measured strains, employing from 8 to 18 gage3 
mounted as sho?'n in figure 10, of refer snce 5. It is likoly 
that the average of these measurements closely represents 
the average strain in the stringers. Hence, the experimen- 
tal results thus obtained should be compared with the cal- 
culated effective width corresponding to the integrated 
mean value of the plate load between nodal lines, rather 
than the load existing at the nodal lines. The effective 
width based on this integrated mean value of the load can 
be o.btained immediately from the integrated mean value of 
equation (30) .that is, 
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fl£) = m ftf 1 E a _ ,1 ■- I 

V> /mean a J lj2(i--u s ) Est F 9 I'd-v 
o -»- 



o 

3 /is 



a )": 



! st 



" fe) ]['io C0S + 'is COsS * x + F 13 8inS * x ]} ftx 

L(i-y 2 )E s t vw ■ J 



2(i-v s ) E at 3P 9 

where I, B = ?ia + gl3 

2 

The value of F , as a function of y., is indicated 

X 5 

in figure 16. Assuming Young's moduli to be the same, for 
the plate and stringer the mean effective width curves for 
the two limiting conditions, that is, p, = 0 and jj, = 00 , 
are shown in figure .17. The.ef f ective width values given 
by equation (32) are slightly lower than those given by 
equation (31) , the maximum difference being of the order 
of 6 percent. 

In calculating the effective width for. any plate- 
stringer combination, it should be kept in mind that the 
torsional rigidity of the stringer varies a B O b % increases. 
S.ince both Pg .. and F 1B are functions .of u, it is neces- 
sary to. correct these quantities for the change in torsion- 
al rigidity of the stringer as O st "increases. Equation 

(25) of reference 5 can be written in the form 

E s t °BT + <°xlp " °> 0 " Sy. Bin Sp= 0 (33) 

where ' ■ • - ■ 

E s t Young 1 s modulus for the stringer 

Ojrp torsion bending constant (references 11 and 12)* 

C x axial compressive stress in the stringer 

t .polar . moment of. inertia of th<e stringer section 
^ about the axis pf twist 

*See appendix' A, for example. 
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0 torsional rigidity of the. stringer 

K (I sin EEL2. moment transferred by the "buckled plate to 
a the stringer 

$ unit torsional deflection of the stringer 

Neglecting end effects, a solution of the above equa- 
tion can he written in the form 

$ = A sin EELE 
from which a 

A = ^-«--r~T7-r-T^-~-7-TT^T 3 (3 4) 



! 8t 0 B ^m!T,JI + (C-a x lJ)j^Lj s 



Since the amplitude A is a measure of the torsional 
stiffness of the stringer the ratio 

Or _ Ao _ E^t O BI (V) + (0 - Ox l-o ) (gB) 
° A E st 0 BT C 

would be a measure of the change of the torsional stiff- 
ness of the stringer, where C r is the reduced torsional 
stiffness of the stringer and A 0 corresponds to A 
when ox = 0. This equation then allows a calculation 
of the variation in torsional rigidity of the stringer 
with increasing stress. 

The experimental values of effective width for a num- 
ber of panels have been compared with the results given 
by equation (32) and are shown in figures 18 and 18a. 
The curves are the calculated values corrected for the 
change in torsional stiffness of the stringe-r with increas- 
ing stress, and taking into account any difference between 
E 8 and E st . As indicated by this comparison, the cal- 
culated results do not compare ' favo rably with the experi- 
mental results. This is particularly true for the torsion- 
ally weak stringer. Although the correction for reduced 
torsional stiffness of the stri-ng-er recognizes a change 
in wave form of the plate with increasing stress, it is 
felt that the magnitude of this correction does not fully 
compensate" for the actual changes in wave form which may 
occur. It is believed that this is the principal cause of 
the deviation of the calculated values. It should also be 
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noted that the derived equations are only valid below the 
proportional limit. It has "been proposed in literature 
on the subject of effective width, that above the propor- 
tional limit it is only necessary to replace E s by E g 

where E s is the secant modulus. This is a fallacy, for 

only that portion of the sheet subjected to stresses above 
the proportional limit is affected by the change in E. 
The distribution of stress in that portion of the sheet 
subjected to stresses below the proportional limit remains 
unchanged. Above the proportional limit it is necessary 
to make a point-by-point correction for the reduced value 
of E. It will be appreciated that a blanket application 
of the secant modulus will reduce the stresses throughout 
the entire sheet and thereby incorrectly reduce the effec- 
tive width. In figures 19 and 20, a point-by-point cor- 
rection is illustrated for the J-section panel with .040- 
24S^-T alclad sheet. The method consist's essentially of a 
calculation of the strain e^' , , at the median fiber in 



the x direction, where 

e„ + ve„ 

■ .(36) 



* 1 - V" 

The working expression for c^' was obtained directly 
from equations (16) and (18) . This strain distribution 
was then plotted and the corresponding stress in the sheet 
obtained from the stress-strain- curves of figures 36. Of 
course, using the stress-strain curve obtained from a uni- 
directional load Is not strictly correct. However, if ay- 
is small in comparison to a x the method should at least" 
approximate the correct stress distribution. A study of 
these curves indicates that the reduced E, above the pro- 
portional limit, is not sufficient to account for the dis- 
crepancies of the effective width calculations indicated 
in figures 18 and 19. Of course, the above method does 
not consider the effects of the induced bending stresses. 
The dimensions of the stringers are shown in figure 42, 
and the stress- strain properties in figures 36 to 40. (The 
data on the J-section was supplied by the Lockheed Aircraft 
Corporation) 

The x components of the bending stress are computed 
as follows: 

EZ M_ 

" d (37) 
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where Z is measured, from the median surface of the plate. 
Substituting for Z its maximum value t/2 and for M x 
the expression 

equation (37) becomes 

°* B_ snorts? + » i?y (38 > 



Substituting the expression (1) for w, equation (38) 
can be written in the form 



tin 0x (39) 



whore 



3T 
is 



|^A(ua 8 - 0 a ) cosh ay - B (vp s + tf 8 ) cos |3yJ 



The function F ie has been evaluated as a function of p, 
for a number of values of" y, as indicated in figure 21. 
If both sides of expression (34) are divided by a cr , it 
can be written in the form 

• Noting that, „ \1/ S E t e 

°cr = ~ — a — 5" — g- 
12(l-v )b 

the above equation can be written as 
2* B _ ■- 1.732 /Eg T 1T / a st ~("i-v S ) 

This ratio assumes its maximum value when sin $x = 1. 

Hence ( -Si J occurs along the line x = and is 

V°cr / max . - . 2m 

given by the bracketed portion of equation (35). 
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The sum of this maximum bending stress distribution 
and tn.e compression stress distribution given by .expres- 
sion (27) gives the maximum compression stress distribu- 

tion in the x direction. Taking = 10 and B s t = 3 S , 

Oct 

the maximum stres-s distribution has been plotted for \x = 0 . 
and p, = oo, a nd is shown in figure 22. As seen from this 
figure, for the case of simply supported edges the maximum 
compression stress does not appreciably exceed the stringer 
stress over the entire panel width. The value of \x = 0 . 5 
corresponds to a stringer of high torsional rigidity and 
in this case, at the edges y = the maximum compression 

stress exceeds the stringer stress by approximately 27 per- 
cent. This increase in the maximum plate compressive 
stress over the stringer stress will not affect to any ap- 
preciable extent the load-carrying ability of the Bheet. 

The wave form of equation (l) was compared vith the 
measured wave forms given in figures 27 and 28 of reference 
4, as shown in figure 23 of this report. Both the trans- 
verse and longitudinal wave -forms check closely with the 
measured ones. In order-to make this comparison, the tor- 
sional rigidity of the stringers must be known. It was 
computed approximately by assuming- the stringer rigidity 
to be the same as that of a rectangular plate having the 
same cross-section area and thickness as the stringer". 
By this method the torsional rigidity was found to' be 630 
inch-pounds per radian per inch. The corresponding value 
of jj, was 0.198, the wave length was determined from fig- 
ure 2 and the maximum amplitude from. equation (26). The 
buckling stress of the sheet, was computed and found to be 
2500 pounds per square inch, which checks with the value 
indicated in figure. 55 of reference 4. 

In view of the good agreement between the calculated 
and measured wave, forms* it is reasonable to expect also 
a good agreement between the calculated -and mea sured 
strains. Calculating th'§ strain by the - equat ion 



e x = 5s. ,+ I'fZsS ■ i- i it* li-ull a /Per] _ ' , 



20: 



cos 2 0x] - -3dL (41) 
E st 
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the resulting values at x = — and ~ are indicated in 

" _ ■ m 4m 

figures 24 and 25, where they are compared with the meas- 
ured values given in figure 20 of reference 4. The meas- 
ured and calculated strain distributions compare favorably. 

The transverse strain at the median fiber, at x = Qjl-L^ ^...a, 

m 

has been calculated and is compared, in figure 25, with 

the observed value at x = fi-i-i. 5 -§._iL. The agreement appears 
to be remarkably good. m 



II- EXPERIMENTAL RESULTS 



Inasmuch as the wave form is one of the important 
factors in the theoretical calculations of the stress dis- 
tribution in a buckled plate, it is desirable to have some 
knowledge regarding its variation with increasing load. 
A number of tests were conducted in which the wave form 
was measured for various increments of load. The measure- 
ments were obtained by means of a dial gage attached to a 
frame mounted on the bed of- the testing machine. The re- 
sults of these tests are shown in figures 26 to 33, where 
the deflections were plotted for various values of the 
ratio of the stringer strain e s t to the sheet strain e 0 
at buckling. The primary purpose of these tests was to 
determine qualitatively the amount of change in the wave 
form with increasing stringer stress. The wave forms plot- 
ted in figure 29 indicate that there is in some cases a 
considerable change in the wave form. The_change in wave 
form as indicated by these curves is somewhat irregular: 
for example, the 0.025 sheet changed but slightly in wave 
form; whereas the change in the 0.04Q sheet iB more than 
in the 0.050 sheet. It seems there should be some degree 
of consistency in the amount of change in the wave form 
as the sheet thickness increases. 

It is possible that the inconsistencies indicated by 
these results can be attributed to the influence of buck- 
les in the adjacent sheet panels. It is difficult to es- 
timate the influence on the compressive stress distribu- 
tion due to the amount of change in wave form as indicated 
by the experiments. However, the indicated flattening of 
the wave will cause an increase in the y component of 
strain which, in turn causes a decrease in the x compo- 
nent of strain. 
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The compression stress-strain properties of the 24S-T 
sheet material used in the above tests are shown in figures 
34 and 35. In addition, a number of compression tests were 
conducted on 24S-T alclad sheet; the curves for this mate- 
rial are shown in figure 36. The method employed in deter- 
mining the compressive stress-strain properties was identi- 
cal to that described in reference 7. 

The stringer strain as a function of average stress 
is shown in figures 37 to 39 for the various panels tested. 
By average stress is meant the compressive load on the pan- 
el divided by its total cro ss- sect ion area. The stress- 
strain properties of the stringers are shown in figures 40 
and 41; the section dimensions are shown in figure 42. 



Guggenheim Aeronautics Laboratory, 

California Institute of Technology, 
Pasadena, Calif., June 1941. 
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APPENDIX A 



The value of the t ors ion- bending constant,. C^y , may 
he calculated either analytically or graphically. Con- 
sidering, a general section such as that shown in figure 
43a, the value of Cbj is given by the expression 

C BT = f w2dA " j[/ w dA ] (42) 
A A 

where w is the circumferential warping and is given by 
the equation 



u 



0 

In the above expression the cross-sectional warping 
r n n has been neglected. Since n, at the most, is 
equal to t/2 , the warping r n n will, in g-ene.ral , be 
small compared to the circumferential warping and may be 
neglected. 

In evaluating the warping w, the following sign 
conventions have been adopted. A direction of rotation 
is chosen arbitrarily; then the circumferential coordinate 
u is taken positive in the direction of rotation. Also 
the positive sense of rotation indicates the posit-ive 
direction of-~the tangent to the 'circumferential coordinate 
u. from the positive tangential direction, the positive 
normal n points to the right and a line drawn from 0 
in the positive n direction indicates the positive r n 
direction. The angle a is measured in the positive 
sense of rotation, that is, from +z to +y. As an illus- 
tration, consider the channel section shown in figure 43b 
and first calculate C-g^. about the point 0. 

Choose the origin of u at 0; then, 

From 1 to 2 r^ = 0 

From 2 to 3 r. = ^ du is negative and w = - — 

t 2 2 

From 1 to 4 r ^ = 0 

From 4 to 5 r^ = & du is positive and w = ~ 

3 3 
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from which 

D 



I* wdA = y udu + Sty* udu = 



6 

o 



Now let Ojrjj "be calculated about the point Qj, 
Again, take the origin of u at 0; then, 



Irom 1 to 2 r t =s + c du is negative, and w = - cu 

ffrom 2 to 3 r t = +f- du is negative, and w = „ 

2 2 



Trom 1 to 4 r^ = +c du is positive, and w "= cu 
From 4 to . 5 r t =. & du is positive, and w = Sfi, + 

<3 '2 



from which 
and a/s 



wdA = 0 

'A 



■ 4 , 

^C. 2 t y u s du + ftjt J (u + e) S duV = %._L ie S ( a + 6l3 ) 

o o v ^ 

+ 2b s (b + 3c) j- 

The integral / wdA will not be zero if, for example, 

A 

the origin of u is taken at point 2. However, the alge- 
braic sum of the two integrals of equation (42) will always 
be constant for a given axis of rotation, .that is,- the ex- 
pression for C 32 j is invariant with respect to the origin 
of the circumferential coordinate u. The origin of the 
rectangular coordinates xy must, of course, be taken "at 
the axis of rotation. Assume, for example, that the cen- 
ter of rotation of the section shown in figure 43c is at 
the point 0; then O^j is evaluated as follows: 

Take the origin of u at point 4; then-from 4 to 3 

r t = r + d cos 8 + b sin 9 . 
and * 6 



*4_3 



■J r t du = r J {t + d cos 8 + b sin 8) d8 



Prom 3 to 2 = r - d 

and 
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u 



w_ = w _ 

3—2 4—3 



+ / ( r — d) du = r 3 Tr + 2rb + (r - d)u 
9=TT </ 



From 2 to 1 r t = 0 

and w . = w j = -r S TT + b(3r - d) 

3_i 3 ~ a j u="b 

Hence ' 

TT 



.A ^ A o a 

+ b] S d8 + t J C*"" + 2rb + (r - d)u] 8 du + t J [r 2 TT 

O ^ TT ° 

+ b(3r - d)] 2 du - I . r s t /[r8 + d sin 8 ~ b cos 8 + b] dR 

b 0 E 

+ t JlT 3 -n + 2rb 4- (r - d)u] du + t [t'tt + b(3r - d) ] du'f 

where A is the total cro sb- sect ion area of the section. 
It can be seen, for example, that, if the origin of the 
circumferential coordinate u is taken at point 1, the 
resulting expression for C-g^ would be Cjons iderably sim- 
plified, since the warping from 1 to 2 vould be zero. 

When a longitudinal stiffener. is attached to a thin 
plate the rigidity of the plate requires that .the a~ie of— 
rotation lie in the plane of the plate.- In computing . 
C-grp for this condition, a number of values of cor- 
responding to various positions in the plane of the plate 
were comput~e~d, and a curve was then plottod of" Cgj as 

a function" of position. The minimum value of was 

then used in comput ing. the variation of torsional rigidity 
with increasing load. 
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Figure 19— Axial stress and strain distribution beyond proportional liuit. h! 



i 



















































































\ 








-1 1 




or const! 


J' ... 




* 






d x corr« 
to 


DBpon 


ding 

V 


/ 






— r- 

\ 














V 

' * J 














/Pr< 


3P. li 


ait./ 
















/ 












/ 






























\ — 

V 








cat 
o/ 














\ in 






A!" 
fly 


















M 


















"v 


/ 


i ■ ■■■ 

1 
1 







W2 "b/2 





... ... 
















/ 

/ 


— \ 


















/ 

/ 












6„ for constfi 
* E 


ait ' 
V 








0' 












— /— 
/ 

/ 


i 












corr( 
to 


jepon 


iing 


i 


















r— H; 
1 

1 


. ,\ _ 


















1 

iiA 




















JTro 


p. lis 


ait^ 








u 








1 












A! 


— _ 





J 

l a. 














V 


— _ 


















1 — \ 


— i 




r 






____ 

































/ 











p 



►3 
o 

I 

O 
P 

*; 
o 

!2i 
O 



CD 

to 



e gT = 70 x 10" 



Ogrn = 50 X 10 

Figure 80.- Axial strain and stress district ion boyond proportional limit 



8 



i 



NACA Technical Sote So. 859 




2TAGA. Technical Hote ITo. 859 



13 



12 



11- 



10 



■>cr 




















^.=0.5 
\ \ 














/ 


\ 
\ 

\ 














/ 
/ 


"T" 

\ 














T 






\ 
\ 

W 






, 

// 




\ 

\ 













// 

y 




\ 

\ 










- / 


7*" 




















- 










i 
















































































i 









b/2 



Figure 22.- Combined axial stress distribution at x= So - 
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Figure 24.- Comparison of calculated axial strain distribution with experimental results given in 
reference 4. 
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Figaro 25.- Comparison of calculated strain distribution with oxporimontal rosults given in rof. 
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Figure 38.- Bull) angle 10S65, 24ST "sheet, paaaol width. = 15 in,, total 
stringer area a 0,1910 sq. in. 
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Figure 39.- Bulb angle 1, 24ST saaet, panol tfidth^ 15 in., total stringer aroa = 0.1130 sq. in. 
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